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Summary 

This is a guide for the use of the pressure disk rotor model that has been placed in the incom- 
pressible Navier-Stokes code INS3D-UP. The pressure disk rotor model approximates a helicop- 
ter rotor or propeller in a time averaged manner and is intended to simulate the effect of a rotor in 
forward flight on the fuselage or the effect of a propeller on other aerodynamic components. The 
model uses a modified actuator disk that allows the pressure jump across the disk to vary with 
radius and azimuth. The cyclic and collective blade pitch angles needed to achieve a specified 
thrust coefficient and zero moment about the hub are predicted. The method has been validated 
with experimentally measured mean induced inflow velocities as well as surface pressures on a 
generic fuselage. Overset grids, sometimes referred to as Chimera grids, are used to simplify the 
grid generation process. The pressure disk model is applied to a cylindrical grid which is embed- 
ded in the grid or grids used for the rest of the configuration. This document will outline the 
development of the method, and present input and results for a sample case. 

Nomenclature 

Note: Unless specified, the following variables are made nondimensional as follows: lengths 
are made nondimensional by radius, areas are made nondimensional by radius , velocities are 

made nondimensional by , densities are made nondimensional by , pressures are made non- 

2 2 
dimensional by , forces are made nondimensional by (U^radius) , and torques are 

2 3 

made nondimensional by p U radius . 
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Aq 

Collective pitch angle, radians 

Ai, 

Cyclic pitch angles, radians 

c 

Rotor blade chord 

Ct 

Thrust Coefficient, — thrust 
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dA 

Elemental area on rotor disk 

dD 

Sectional drag of rotor blade 

dL 

Sectional lift of rotor blade 
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Sectional torque of rotor blade 

Elemental span 

Sectional thrust of rotor blade 

Azimuthal spacing of grid, radians 

Number of blades 

Unit vector normal to blade path 

Unit vector parallel to blade path 

Pressure jump 
Rotor radius 
Radius location 
Radius vector 
Radius unit vector 
Freestream velocity, ft/sec 

Velocity magnitude relative to rotor blade 
Velocity vector 

Induced component of velocity normal to the rotor disk 
Tangential velocity jump 
Rotor blade tip velocity 

Azimuthal component of velocity relative to rotor blade 

Normal component of velocity relative to rotor blade 

Effective angle of attack, radians 

Geometric angle of attack, radians 

Induced angle of attack, radians 

Geometric blade twist, radians 

V- 

Mean induced inflow ratio, ^ 

tip 

Advance ratio, ^ 

tip 

Focal air density 

Free stream air density, slugs/ft 

Azimuth, radians 

Shaft rotation vector, radians/sec 


Introduction 

The flow field in which rotorcraft operate is dominated by vortical flow in the rotor wake and 
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viscous effects at the fuselage. Early analytical work focused on development of models to pre- 
dict the performance of the isolated rotor without consideration of the effect of the fuselage. 
However, the fuselage can have a significant effect on the overall performance of the rotorcraft. 

The goal of this code is to model the effect of the rotor on the fuselage and also give an indi- 
cation of the flow field at the rotor. By accepting a time averaged solution for the rotor and wake, 
the complexity of the problem is reduced, while still yielding a time averaged rotor influence on 
the fuselage. With this method it is not necessary to generate grids for the individual rotor blades, 
reducing the cost in terms of operation and grid development and allowing quick changes in the 
rotor geometry to be made without generating new grids. This method will also provide the user 
with an estimate of the flow at the rotor and the rotor load distribution. 

The method presented here uses overset, or Chimera, grids (Ref. 1) to allow the grids for the 
rotor disk and the fuselage to be generated independently. The incompressible Navier-Stokes 
code INS3D-UP (Ref. 2) is used to provide a flow solution. While the flow field relative to the 
rotor blade motion is certainly compressible, it is assumed that the primary effect of this is in the 
magnitude of the forces generated by the rotor blades. Thus, the effects of compressibility can be 
included when determining the forces generated by the rotor blades. 

Pressure Disk Model 

The method presented in this paper simulates the helicopter rotor by imposing boundary con- 
ditions at the surface swept by the rotor in space. The boundary conditions are applied on adja- 
cent planes of a cylindrical grid. At this surface, a disk, two separate boundary conditions are 
imposed. The top of the rotor disk surface is modelled with an outflow boundary condition. On 
this plane, the pressure and velocities tangential to the disk are extrapolated from the flow field 
above the rotor. The velocity normal to the disk is specified to be equal to the normal velocity at 
the bottom of the disk to maintain continuity. The lower disk surface is modelled with an inflow 
boundary condition. On this plane, the pressure and tangential velocities are based on their values 
at the top of the disk and the normal component of velocity is extrapolated from below the disk. 
The effect of the rotor is imposed by the difference in pressure and tangential velocities between 
the upper and lower disk surfaces. 

The first step in calculating the rotor model boundary condition is to determine the velocity of 
the flow relative to a rotor blade laying in the rotor disk plane. Once the flow at the disk is known, 
it is possible to determine the forces generated by a rotor blade at each point on the disk. Blade 
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element theory is used to compute the rotor forces. 

The flow solver provides the velocity, V , relative to the fixed frame axis at every point on the 
disk grid. Using Fig. 1 as a guide, it is possible to determine the velocity at point 5 on a rotor 
blade whose path coincides with the disk grid. Using the grid definition, the unit vector normal to 
the blade path, n , can be determined. V ■ n is the velocity of the air in the direction normal to the 
blade path, f is the vector from the shaft center to the point B . Because r is perpendicular to n , 
the vector parallel to the path of point B, p , can be determined from r and n . 

_ ^ ^ f' 

p = rxn; f = rr, (1) 

\r\ 

V ■ p is the component of the freestream velocity in the direction parallel to the blade path at 
point B. (r X Q.) ■ p is the velocity due to the blade rotation in the direction parallel to the blade 
path at point B . The total velocity parallel to the blade path then becomes V ■ p + (f-xH) ■ p. 

Using the velocity components normal and parallel to the blade path, the induced angle of 
attack, a . , can be determined. 


f V-n \ 

a. = atan :: 

(2) 

' ■ p + (rx Q.) ■ p'' 

The effective angle of attack, a^, at point B is 



(3) 


The relationship between a , a , and a. is shown in Fig. 2. a is the geometric angle of 
attack and is determined by blade twist and the pitch of the blade. The pitch of the rotor blade rel- 
ative to the shaft of the rotor is given by 

= AQ-A^cos\|/-BjSin\|/ -I- 1^0.75 . (4) 

Once the velocity relative to a blade element has been determined, the lift and drag generated 
by the blade element can be calculated. The grid used for the rotor is cylindrical and a planform 
view is shown in Fig. 3. Rotor geometry will specify that the rotor blade passing through the area 
dA has a chord of c. The grid definition will determine the value of the elemental span, ox dr . 
With the effective angle of attack calculated, values for the lift and drag coefficients can be deter- 
mined. These values could come from airfoil tables, other computational methods, or other sim- 
ple empirical models. The user is given two options for calculating the blade forces. The first 
option is to set the lift coefficient to With this option, once the lift coefficient reaches a 

user specified limit, it is not allowed to increase. While this is a rather simple stall model, it has 
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proven reasonably accurate. The drag coefficient is also determined with a simple empirical 

2 

model defined as 0.006 + 0.01 33 . In the second option, the model uses airfoil tables to deter- 
mine the forces. Because of the added expense of searching and interpolating the data in the air- 
foil tables, this option is slightly more expensive to use. If airfoil tables are to be used, they must 
be in the C81 format (Ref. 3). The variable model in the input for the rotor boundary condition 
controls this option. The amount of lift and drag acting on the section of the rotor blade passing 
over the area dA is 


dL = c,cdr 
dD = . 


(5) 

( 6 ) 


In reality, the area, dA , would not see the effects of the rotor blade continuously. For this rea- 
son, the added pressure and velocity jumps are scaled by a time factor, where N is the 
number of rotor blades and dtj/ is the grid azimuth resolution. This is the traction of time a rotor 


blade spends over the area dA . With this time factor, dL, dD, and the forces can be resolved 


into the thrust and torque components. 

dt|/ 


dT 


(dLcos (a ) - dDsin (a ) ) 
271 ' ' 


dQ = (dLsin (a ) -i- dDcos (a ) ) 
2tc ' ' 

The pressure jump boundary condition then becomes 


(7) 

( 8 ) 


A d.T 

AP = ^ . (9) 

The torque boundary condition is modelled as a jump in the tangential velocity. The magni- 
tude of this jump is based on the induced velocity, and the ratio of torque to thrust. The 
induced velocity was calculated when the velocity relative to the rotor blade was determined. 

Av = V. (10) 

In practice, the torque boundary condition effect has been smaller than expected and in some 
cases its use has led to nonphysical solutions. Other, potentially more effective, ways of imposing 
the torque boundary condition are being studied. It is not implemented in the current version of 
the model, but it is discussed here for completeness. 


Trim Loop 

At user specified iterations, the model calculates the correct collective and cyclic pitch angles 
in order to achieve the desired thrust coefficient and eliminate moments about the hub. In the cal- 
culation of the moments, the flapping hinge offset is assumed to be zero. The trimming routine 
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assumes linear relationships between the thrust coefficient and the collective pitch angle, and 
between the hub moments and the cyclic pitch angles. A simple linear relationship is assumed to 
exist between the control angle and the error that is being minimized. Figure 4 shows the general 
scheme used to implement the trim loop. The first step is to calculate the pressure jump used for 
the boundary conditions. This can be integrated over the disk surface to determine the thrust 
coefficient and the flapping moments about the hub. Once this is done, the trim subroutines 
search for the proper pitch angles. The pressure jump is then recalculated using these new angles. 
This loop continues until the new angles match the angles from the previous loop, or a user speci- 
fied number of iterations has passed. The Navier-Stokes solution is then advanced with the new 
angles. When several trim cycles pass with no change in the control angles, the rotor model is 
assumed to be converged. This does not necessarily mean that the overall flow field solution has 
converged. 

Appendix A discusses the variables used by the rotor model boundary conditions contained in 
the file " rotor inp\ Several of these variables are used to control the trimming operation. Initially, 
the collective and cyclic angles can be set by the user. If the user does not wish the model to trim 
the rotor, itrmst should be made sufficiently large so that the trimming routines are not called. 
It is not possible to change the collective and cyclic pitch angles at a restart. The only means of 
adjusting these angles on a restart is through the trim routines. The variable iramp is used to 
specify the number of iterations to ramp up the pressure jump. This can be used if the required 
pressure jump is so large it causes difficulty with the convergence of the code. The variable 
iclcf rc controls the intervals at which the rotor forces (C-p and hub moments) are calculated 
and printed. This is independent of the trimming operation, but does provide an insight to the sta- 
tus of the solution. The variable inttrm controls how often the trim subroutines are called. For 
difficult cases, inttrm should not be too small as it is best to allow the solution to “calm down” 
before trimming again. The variables altol, bltol, and cttol are the tolerances used for the 
trimming. If the trim routines determine the change in the hub moments and thrust coefficient are 
smaller than these values, the trim routines are exited. Setting altol and bltol to ‘-888’ will 
prevent the cyclic pitch angles from being adjusted when the boundary condition is used to model 
a propeller. The variables aline, blinc and ctinc, in degrees, are used to keep the trim- 
ming routines from diverging during the early attempts to trim by setting the maximum change in 
the control angles. 
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Flapping 

The flapping motion of the rotor is not calculated by the code. However, if the flapping 
motion is known the code can account for the first two harmonics of flapping by adding the flap- 
ping velocity to the velocity normal to the blade path. Coning is handled by using a predeter- 
mined coned grid geometry. 

Computational Method 

As stated earlier, the pressure disk rotor model is implemented in version 1.8-19 of INS3D- 
UP. While this document is not intended to be a manual for INS3D, there are some specific steps 
that can be taken to increase the chances of success in the operation of the code with the rotor 
model. A sample input for the isolated rotor case is included in Appendix B. The most important 
variables to the operation of INS3D with the rotor model are n j swp, nkswp, and nlswp in 
namelist zonin. nkswp should be set to zero, while typically, n jswp and nlswp are set to a 
value of 5-10. A higher number of sweeps can aid in the convergence, but the cost of running the 
code is directly related to the number of sweeps performed. It is also important that kpr is set to 
‘T’. This variable controls the periodic boundary condition in the k direction. Finally, it is not 
possible to change the angle of attack during a restart. 


Grids 

The pressure disk model uses a cylindrical grid, an example of which is shown in Figure 5. 
The generation of this grid can be done in a relatively simple code. Using the jkl convention, the 
plane of the disk should lie on a constant j-plane. The k=l plane is located on the 0° azimuth. 
Uniform spacing must be used for the distribution in the k direction. The k=k m^y plane is located 


at the 


fk - n 
max 


360° azimuth. The 1 direction is such that 1=1 corresponds to the center of the 

V ‘^max J 

grid and l=lniax the outer cylinder of the grid. Note that the 1=1 plane is a singularity. This 
requires the use of the axis boundary condition. The rotor model boundary conditions are written 
so that positive thrust will result in thrust in the -j direction. In other words, the j value at the top 
of the disk is less than the j value at the bottom of the disk. 

As stated earlier, overset grids are used to simplify the grid generation process. With overset, 
or Chimera, grids, complex geometries are broken down into simpler component geometries. 
Individual grids can then be generated for these simpler shapes. Then the collection of grids can 
be combined to produce the final overset grid system. 
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Figure 6 shows an overset grid system for an isolated rotor case. Two grids were used to form 
this system. The first grid is the cylindrical grid used for the rotor model. The second grid is a 
simple cartesian grid that serves to move the far field boundaries sufficiently far from the rotor 
disk. The cartesian grid also makes applying the far field boundary conditions easier. PEGSUS 
4.0 (Ref. 4) was used to combine the individual grids into the single overset grid system. It is 
assumed that the user of this code has experience with PEGSUS. 

A problem arises when using PEGSUS with the cylindrical grid used for the rotor model. Due 
to the manner in which INS3D-UP treats periodic boundary conditions, the cylindrical grid is not 
completely closed, i.e., a slice is removed. When PEGSUS performs a search to determine where 
to cut holes in the outer cartesian grid, it will fail to cut holes at the missing slice. It is important 
to remember that the gap in the cylindrical grid is not seen by INS3D-UP. Therefore, something 
must be done to cut holes in the cartesian grid. This problem is solved by adding an additional k 
plane to complete the cylindrical grid and cut holes in the cartesian grid. In other words, on this 
new cylindrical grid, the k=l and the k=kmax lines coincided. PEGSUS is then run using this 
new cylindrical grid. The file containing the interpolation coefficients is then searched to find 
instances where the kmax grid index is used. In these cases, the kmax value is replaced with 1. 
The new cylindrical grid is then replaced with the old cylindrical grid when the final grid file is 
generated. 

The cylindrical rotor disk grid should be generated with the experimentally determined, or 
estimated value of coning. During the validation of the method, although several cases requiring 
slightly different coning angles were used, only one coning angle was used. Experience has 
shown that the effect of coning on the solution is small and having a grid that is 1 degree off in 
coning makes no noticeable difference in the calculation. 

The sample isolated rotor grid system uses a cylindrical rotor disk grid with 58 points in the 
vertical direction, 50 points in the circumferential direction, and 68 radial points. The outer carte- 
sian grid has 65 points in the streamwise direction, 65 points in the lateral direction and 50 points 
in the vertical direction. Spacing in the vertical direction is approximately 0.0003 R at the disk 
surface. Geometric spacing is used above and below the rotor in the vertical direction. A cosine 
distribution is used for the radial spacing on the disk with geometric spacing used off the disk. 
The total number of points was 408,450 and 16 million words of memory were required for the 
sample isolated rotor computation. 
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Implementation of Boundary Conditions 

The boundary conditions added to INS3D-UP by the pressure disk model are 101 and 102. 
Boundary condition 101 is used for the top of the disk, 102 is used for the bottom of the disk. The 
top and bottom planes of the disk are separated by a plane with a boundary condition of - 1 
(iblanked). This is an INS3D-UP requirement. Because the first plane in the 1 direction of the 
cylindrical grid is singular, the axis boundary condition (71) must be used on this plane. The 
boundary condition file bcmain.dat for the isolated rotor test case is shown below. The rotor 
boundary conditions must be listed last in bcmain.dat. This does not mean that the grid used for 
the rotor boundary conditions has to be the last grid in the system. INS3D allows the boundary 
conditions to be listed in any order. Appendix C shows a sample of the input file rotor.inp used by 
the rotor model boundary condition. 




c 

c- 

c 

c- 


boundary conditions for isolated rotor 
ibcval zone jbeg jend kbeg kend 


Ibeg lend 


c cartesian 

grid boundary 

conditions 





31 

1 

65 

65 

1 

65 

1 

50 

c rotor grid 

boundary 

conditions 





71 

2 

1 

58 

1 

50 

1 

1 

101 

2 

28 

28 

1 

50 

9 

50 

-1 

2 

29 

29 

1 

50 

9 

50 

102 

2 

30 

30 

1 

50 

9 

50 




Sample Results 

The pressure disk method was tested with the rotor alone as well as with the addition of a 
generic fuselage at five advance ratios, |i, ranging from 0.15 to 0.40. Results of these tests were 
compared with experimental induced inflow measurements and trimmed control angles (Refs. 5- 
9). A more extensive discussion of these results can be found in Ref. 10. The trimmed control 
angles for the isolated rotor test case with both airfoil models are compared with the experimental 
values for the 0.15 advance ratio case reported by Ghee and Elliot (Ref. 11) in Table 1. 


Table 1 : Computed Trim Conditions 


Advance 
Ratio, p. 

Experiment 

Isolated 

Rotor 

w/lTta 

model 

Isolated 

Rotor 

w/airfoil 

tables 

0.15 

Ao=6.55° 

Aq=6.76° 

Ao=5.76° 


A;=-1.39° 

Aj=-2.2T 

A;=-2.24° 


Bi=\.99° 

B,=2.A2° 

B;=1.80° 
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Code Output 

Two types of information about the rotor is included in the standard output from INS3D-UP. 
The calculated thrust coefficients and the lateral and longitudinal hub moments are printed every 
iclcf rc iterations. This information is useful for determining the rate of convergence of the 
trimming routines. Typically, iclcf rc is set equal to inttrm. The other information also 
relates to the trim routines. At every inttrm iterations, the code prints the number of iterations 
needed by the trim routines, as well as the values for Aq, Aj, and Bj. At the start of the calcula- 
tion, it is not unusual for the number of trim iterations to be equal to itloop. As the thrust coef- 
ficient approaches the value specified in cthrust, the number of iterations should decrease, 
eventually becoming 0. A sample of the standard output in included in Appendix D. 

The code will also give information at the disk. This information is contained in the file 

fort.35. This file is printed every np3d iterations and at the end of each run. The file is formatted 

for use with the Tecplot®(Ref. 12) plotting software, but the format is easily translatable for use 

with other plotting software. The file contains the following information for each grid point on 

the disk: the x and y coordinates, the sectional lift and drag coefficients, angle of attack, induced 

angle of attack, mean induced inflow ratio at the j= jinf Iw plane, the pressure jump across the 

disk and the two components of velocity. The velocities, and , are made nondimensional 

by the freestream velocity. The dimensional sectional lift can be calculated as 

] f 2 2"\ 2 

^ ^ \chord ■ 0 ■ U , (11) 

2 /V t|/ n y oo ’ ^ ^ 

where is the lift coefficient, chord is the physical chord of the rotor blade, is the free 
stream density, and is the free stream velocity. 

Computational Expense 

The number of iterations required for the rotor model to converge depends on the conditions 
which are being modelled. The basic scheme for using the rotor model is to start with a small 
amount of collective and no cyclic controls. After some time, about 10-30 iterations depending 
on the severity of the case, the trim routines are called and the collective will start to increase to 
achieve the desired thrust coefficient. The trim routines will also adjust the cyclic controls to try 
to maintain zero moments at the hub. The amount of change during a single trim cycle is limited 
for stability reasons. Experience shows that this method arrives at the desired conditions more 
quickly than starting at an initial guess of the final conditions. This is partly due to the inability of 


10 



the code to handle abrupt changes in the rotor boundary conditions from one iteration to the next. 

Higher advance ratios typically converge fairly quickly. This is caused by the increase in average 

pressure jump required for lower advance ratio cases. The average pressure jump is approxi- 
2 

mately equal to C/\i , therefore as the advance ratio decreases, the average pressure jump 
increases, and the rotor boundary condition has a greater impact on the flow field and in most 
cases slows the convergence. Table 2 shows the number of iterations required for a typical iso- 
lated rotor calculation. The speed of the code on a Cray C-90 was approximately 160x10 ^ sec- 
onds/iteration/grid point. 


Table 2: Convergence for isolated rotor case 


Advance 
Ratio, jO. 

Iterations 

Cray-90 hours 

Memory 

(Mwds) 

0.15 

180 

3.25 

16 
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Figure 1. Rotor blade. 



Figure 2. Rotor section with relative velocities and forces. 
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Normal INS3D-UP Operation 


Calculate Boundary Condition; 


Advance Flow Solution 



Rotor Model Boundary Condition 


Calculate Blade Forces, 
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Moments about Hub, 
Rotor Boundary 
Conditions 


Adjust Collective 
and Cyclic Pitch “ 
Angles 


Figure 4. Implementation of rotor model and trim loop. 
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Figure 5. Cylindrical rotor grid topology. 


Figure 6. Isolated rotor embedded grid, 
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Appendix A 

Input for INS 3D Rotor Model 

The following variables are used to describe the geometry of the rotor as well as the condi- 
tions for which the rotor is being run. The input is in namelist rtrstf in the file rotorinp. 
Default values are shown in [] . 


rmu advance ratio [0.4] 

aao collective (Aq) (degrees) [0.1] 

aaal cyclic (degrees) [0.0] 

abl Bj cyclic (degrees) [0.0] 

aao, aaal, and abl can not be changed at restarts 

note convention: theta = aao - aaal*cos(t|/) - abl*sin(\|/) 

f al first harmonic cosine flapping (degrees) [0.0] 
fbl first harmonic sine flapping (degrees) [0.0] 
f a2 second harmonic cosine flapping (degrees) [0.0] 

fb2 second harmonic sine flapping (degrees) [0.0] 

note convention: flapping = coning + f al*cos(t|/) + fbl*sin(\|/) + 


btwst 

aspr 

bldnum 

cdmax 

clmax 

clmin 

iramp 

iclcf rc 

itrmst 

inttrm 

altol 

aline 

bltol 

blinc 

cttol 

ct inc 

it loop 


amount of linear twist in blade (degrees) [-8.0] 
aspect ratio of blade, (blade radius/blade chord) [13.0] 
number of blades [4.0] 

upper limit allowed for c^j (used when model = 1) [0.3] 

upper limit allowed for Cj (used when model = 1) [1.2] 

lower limit allowed for Cj (used when model = 1) [-1.2] 

ramp up pressure jump for first iramp iterations [10] 

interval to calculate and print rotor forces [5] 

iteration to start trim calculations [10] 

interval for trim calculations [5] 

tolerance for longitudinal trim [0.0001] 

maximum change in A^ during trimming (degrees) [0.1] 

tolerance for lateral trim [0.0001] 

maximum change in during trimming (degrees) [0.1] 

tolerance for thrust coefficient trim [0.000001] 

maximum change in Aq during trimming (degrees) [0.1] 

maximum number of trimming iterations [10] 
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cthrust desired thrust coefficient [0.0064] 

j inf Iw j plane at which the mean induced inflow ratio is to be calculated [15] 
model model used for sectional loads, 1 = 27ta, 2 = airfoil tables [1] 


Notes: 

1. Typically, iramp is set to 10 to give the flow a chance to develop before trimming starts. 
Use an initial guess of 0.1 for aao and 0.0 for aaal and abl. Doing this and using the 
default limits on the trim loops allows the solution to gradually develop and usually results 
in quicker convergence than making initial guesses for aao, aaal, and abl. 

2. Limits on cl are to approximate stall and give more reasonable results. 

3. Make sure itrmst is greater than iramp. 

4. If trim loops have trouble converging, decrease aline, blinc, and ctinc and/or 
increase inttrm. 
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Appendix B 

An input file to INS3D for the isolated rotor case is show below. Note that iover is set to 
“1” indicating that an overset grid is being used. Also notice that the namelist zonin is repeated 
for each grid. The angle of attack is set in namelist geomin. It is not possible to change the 
angle of attack during a restart. 


$datain 
ntmax=500, 
iover=l , 
iss = l , 
istart=0 , 
iturb=2 , 
ivis=3 , 
beta=l 0 0.0, 
dtau=l . e24 , 
dt=l . e24 , 
reynum=2 . 035e6, 
$end 
$ zonein 
kpr= f, 
iaxsym=l 1 1 , 
n jswp=0, 
nkswp=0 , 
nlswp=5 , 
cdiss=l . 0 , 

$end 
$ zonein 
kpr= t, 
iaxsym=l 1 1 , 
n j swp=5 , 
nkswp=0 , 
nlswp=5 , 
cdiss=l . 0 , 

$end 
$geomin 
alpha=-3 .00, 
$end 
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Appendix C 

Below is the input for a typical run of the rotor model boundary conditions. This information 
is contained in the file rotorinp. 

$rtrst f 
rmu = 0.15, 
aspr = 13.0, 
cthrust = 0.0064 
iramp = 1, 
itrmst = 5, 
aao = 0.1, 
aaal = 0.0, 

abl = 0.0, 

altol = 0.0001, aline = 0.1, 
bltol = 0.0001, blinc = 0.1, 
cttol = 0.00001, ctinc = 0.1, 
inttrm = 5, 
model = 2, 
jinflw = 12, 

$end 
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Appendix D 

Below is output from a typical run of INS3D-UP using the rotor model boundary conditions. 
The parts of the output that are due to the use of the rotor model are shown in bold font. Com- 
ments are shown in { } . 

INS3D-UP Version 1.80-19 
06/06/95 13:45:45 on sabre 

■k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 


Grid Information 
Opening file named: xyz.dat 

This is an unformatted, multiple zone grid file with an iblank array. 

nz = 1 jmax = 65 kmax = 65 Imax = 50 

nz = 2 jmax = 58 kmax = 50 Imax = 68 

total number of grid points = 408450 

largest single dimension = 68 


Rotor input 

&RTRSTF RMU = 0.15, AAO = 0.1, AAAI =0., ABl = 0., 
FA2 = 0 . , FB2 = 0 . , 

BTWST = -8., ASPR =13., BLDNUM =4., CDMAX = 0.3, 
- 1 . 2 , 


FAl = 0 . , FBI = 0 . , 
CLMAX = 1.2, CLMIN = 


IRAMP = 1, ICLCFRC = 5, ITRMST = 5, INTTRM = 5, AITOL = l.OE-4, BITOL = 
l.OE-4, 

AlINC =0.1, BlINC =0.1, CTTOL = l.OE-5, CTINC =0.1, 

ITLOOP = 10, CTHRUST = 6.4E-3, JINFLW = 12, MODEL = 2, SEND 

{The above input shows that the minimum ramping is being used, the collective is starting out 
small, and the model will start trimming on the fifth iteration. } 


&DATAIN NITER = 1, NP3D = 100, NQFILE = 100, NSUB = 1, NTMAX = 50, IGRID 
= 0, IINTRP = 0, INGRID = 1000, lOVER = 1, 

IPATCH = 0, ISS = 1, ISTART = 0, ISUBIO = 1, ITURB = 2, IVIS = 3, BETA 
= 100 . , DTAU = 1 . OE+24, 


DT = 1 . 

OE+24, 

EIGEPS =0.1, EPSCON = l.OE-4, 

PIN 

= 1.5, 

POUT 

= 1., 

REYNUM 

= 2035000 . , 

SEND 





&ZONEIN 

JMAX 

= 3, 

KMAX = 3, LMAX = 3, KPR = 

F, 

lAXSYM = 

111, 

IFLXODR 

3, NJSWP 

= 0, 

NKSWP 

= 0, NLSWP = 5, NTJSWP = 

0, 




NTKSWP 

= 0, 

NTLSWP 

= 2, CDISS =1., SEND 





&ZONEIN 

JMAX 

= 3, 

KMAX = 3, LMAX = 3, KPR = 

T, 

lAXSYM = 

111, 

IFLXODR 

3, NJSWP 

= 5, 

NKSWP 

= 0, NLSWP = 5, NTJSWP = 

0, 




NTKSWP 

= 0, 

NTLSWP 

= 2, CDISS =1., SEND 





Airfoil 

table 

read from file naca0012.c81 





Allocated 

14.296 

million words of memory. 





&GEOMIN 

ALPHA = -3. 

, REFAREA = 1 . , SEND 






Boundary Conditions: 

Number of BC regions read in from bcmain.dat = 5 


Outflow boundary 
50 

axial singularity 
1 

Actuator disk inflow 

ibcval = 

31 

nz = 

1 

j = 

65 

65 

k= 

1 

65 

1= 

ibcval = 

71 

nz = 

2 

j = 

1 

58 

k= 

1 

50 

1= 

ibcval = 

101 

nz= 

2 

j= 

28 

28 

k= 

1 

50 

1= 
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9 50 


Hole region 
50 

ibcval = -1 nz= 

2 j= 

29 

29 

k= 

1 

50 1= 

Actuator disk 

outflo ibcval = 102 nz= 

2 j= 

30 

30 

k= 

1 

50 1= 


9 50 

{NOTE: the rotor model boundary conditions must be in the above order and must be the last 
boundary conditions specified. } 

Total of 0 no-slip walls processed. 

Allocated 111680 words for bcozne arrays. 


Read in 

overlap 

zonal 

info, nz. 

ibpnts : 

1 9792 

Read in 

overlap 

zonal 

info, nz. 

ibpnts : 

2 1376 


Test 

zone 

of 

1 

overlaid 

zonal 

interface 

reports 

a max error of 

0.3464E-13 

for 

base 

Test 

of 

overlaid 

zonal 

interface 

reports 

a max error of 

0. 6566E-11 

for 

base 


zone 2 

(Make sure the above numbers shown in italics are small. They represent the errors caused by the 
interpolation of grid information. This is easy to overlook, but very important. } 

Zone: 1 Max Jacobian: 0.14972E+04 Min Jacobian: 0.12762E+02 

Zone: 2 Max Jacobian: 0.36905E+08 Min Jacobian: 0.17840E+03 

Allocated 1.634 Mwords for turbulence model. 

Baldwin-Barth Tubulence Model 
Version 1.27 


karman constant = 0.41 

emu 

= 9.E-2 

cle 

= 1.2 

c2e 

= 2. 

sigma e 

= 0.70041 

aplusl 

= 26. 

aplus2 

= 10. 

dt model 

= l.E+24 

re_t inf 

= 1. 

Transition 

locations : 


nz, jtranlo, jtranup = 
nz, jtranlo, jtranup = 
Line-relaxation sweeps: 
nz , nt j sp, ntksp, ntlsp 
nz , nt j sp, ntksp, ntlsp 
Walls : 


6666667 


1 , 2*0 

2 , 2*0 

= 1 , 2 * 0 , 

= 2 , 2 * 0 , 


2 

2 


Distance function was read from file distfun.dat. 


resmaxO for zone 1 
resmaxO for zone 2 


6.76634930277088E-22 
9 . 876532829657071E-9 
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nt resmax divmax j k 

cdp clp 

1 

nz 

turres cd 

cl 








1 0.988E-08 0.988E-10 31 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

27 

50 

2 

0 . 883E-17 

0 . OOOE+00 

2 0.353E+04 0.353E+02 30 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

46 

51 

2 

0 . 916E-04 

0 . OOOE+00 

3 0.322E+04 0.322E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

29 

51 

2 

0 . lOlE-03 

0 . OOOE+00 

4 0.240E+04 0.240E+02 32 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

48 

46 

2 

0 . 931E-04 

0 . OOOE+00 


-7.562800129194934E-4, -1 . 242615678921699E-2 , 4 . 478379790770992E-2 

right left front back 

-0.272E-01 -0.396E-01 -0.176E-01 -0.624E-01 

starting trim loop 
starting trim loop 
10 trim iterations 

AO =1.099999999999987 A1 =-0.3861963347996582 B1 =0.2735704129524734 

5 0.227E+04 0.227E+02 31 27 50 2 O.lllE-03 O.OOOE+00 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

{After the 5th iteration, the calculated Cj is -T.SdxlO'"^, the trim loops are started and after 10 
passes through the trim loops, the collective has been increased to try to correct for the low Cj. 
Note that the cyclic angles have also changed. Had CTINC not been set to a small number, Aq 
might have become too large during the early trimming attempts and the solution might have 
become unstable. } 


6 0.265E+04 0.265E+02 31 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

29 

50 

2 

0 . 147E-03 

0 . OOOE+00 

7 0.288E+04 0.288E+02 31 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

29 

50 

2 

0 . 194E-03 

0 . OOOE+00 

8 0.308E+04 0.308E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

48 

51 

2 

0 .216E-03 

0 . OOOE+00 

9 0.250E+04 0.250E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

46 

51 

2 

0 .209E-03 

0 . OOOE+00 


5.230469709086374E-4, -5 . 173697433143054E-3, 4 . 058687632199098E-2 

right left front back 

0.144E-01 0.922E-02 0.239E-01 -0.167E-01 

10 trim iterations 

AO =1.989999999999938 A1 =-0.7501364027702522 B1 =0.2835252502170498 


10 0.185E+04 0.185E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

44 

51 

2 

0 .230E-03 

0 . OOOE+00 

11 0.155E+04 0.155E+02 31 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

28 

50 

2 

0 .249E-03 

0 . OOOE+00 

12 0.187E+04 0.187E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

29 

51 

2 

0 . 444E-03 

0 . OOOE+00 

13 0.218E+04 0.218E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

29 

51 

2 

0.553E-03 

0 . OOOE+00 

14 0.191E+04 0.191E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

31 

51 

2 

0.715E-03 

0 . OOOE+00 

1.508858882341035E-3, -1 . 588902451108298E-2 , 

right left front back 

3.13936377128674E-2 


0.506E-01 0.347E-01 0.523E-01 0.209E-01 

10 trim iterations 

AO =2.879999999999882 A1 =-1.027095750198953 B1 =0.3014633315327107 


15 0.207E+04 0.207E+02 31 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

29 

50 

2 

0 . 695E-03 

0 . OOOE+00 

16 0.196E+04 0.196E+02 31 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

29 

50 

2 

0 . 647E-03 

0 . OOOE+00 

17 0.234E+04 0.234E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

33 

51 

2 

0 . 122E-02 

0 . OOOE+00 
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18 0.267E+04 0.267E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

35 

51 

2 

0.763E-03 

0 . OOOE+00 

19 0.250E+04 0.250E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

35 

51 

2 

0 . 656E-03 

0 . OOOE+00 


2 . 595043126406765E-3, -2 . 437389306736159E-2 , 1 . 943479592517771E-2 

right left front back 

0.864E-01 0.620E-01 0.864E-01 0.669E-01 

10 trim iterations 

AO =3.76999999999984 A1 =-1.199683126443517 B1 =0.6282236783746846 


20 0.232E+04 0.232E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

35 

51 

2 

0.510E-03 

0 . OOOE+00 

21 0.215E+04 0.215E+02 31 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

33 

50 

2 

0.367E-03 

0 . OOOE+00 

22 0.336E+04 0.336E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

37 

51 

2 

0.459E-03 

0 . OOOE+00 

23 0.332E+04 0.332E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

37 

51 

2 

0.345E-03 

0 . OOOE+00 

24 0.279E+04 0.279E+02 31 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

35 

50 

2 

0.393E-03 

0 . OOOE+00 

3.658183123326649E-3, -1 . 495913391743819E-2 , 

right left front back 

1 . 626192560556116E-2 


0.113E+00 0.979E-01 0.125E+00 0.109E+00 

10 trim iterations 

AO =4.659999999999798 A1 =-1.357818367685731 B1 =0.6576795697164393 


25 0.282E+04 0.282E+02 31 35 50 2 0.418E-03 O.OOOE+00 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

nt resmax divmax j k 1 nz turres cd cl 

cdp clp 


26 0.274E+04 0.274E+02 31 37 50 2 0.447E-03 O.OOOE+00 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

27 0.255E+04 0.255E+02 31 37 50 2 0.651E-03 O.OOOE+00 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

28 0.196E+04 0.196E+02 31 37 50 2 0.663E-03 O.OOOE+00 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

29 0.171E+04 0.171E+02 26 12 38 2 0.701E-03 O.OOOE+00 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

4.817764003763158E-3, -3 . 107120529424989E-2 , 2 . 069634084709637E-2 

right left front back 

0.157E+00 0.126E+00 0.170E+00 0.150E+00 

10 trim iterations 

AO =5.45100106978208 A1 =-1.560228468702284 B1 =1.128201772374972 

30 0.183E+04 0.183E+02 26 12 38 2 0.803E-03 O.OOOE+00 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

{After 30 iterations, the trim routines are still requiring the preset limit of 10 iterations. This 
means that they are not able to trim the rotor. However, note that the C-p is now 0.00482, the 
moments about the hub are under control, and the residual (resmax) is generally decreasing.} 


31 0.179E+04 0.179E+02 26 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

12 

38 

2 

0 . 967E-03 

0 . OOOE+00 

32 0.201E+04 0.201E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

39 

51 

2 

0 . 115E-02 

0 . OOOE+00 

33 0.188E+04 0.188E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

40 

51 

2 

0 . 132E-02 

0 . OOOE+00 

34 0.226E+04 0.226E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

37 

51 

2 

0 . 129E-02 

0 . OOOE+00 


5.753871723164522E-3, -1 . 513962421256565E-2 , 2 . 29005854214348E-2 

right left front back 

0.178E+00 0.163E+00 0.208E+00 0.185E+00 

10 trim iterations 

AO =5.799686119039791 A1 =-1.765531638217212 B1 =1.372259070130418 
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35 0.361E+04 0.361E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

37 

51 

2 

0 . 124E-02 

0 . OOOE+00 

36 0.451E+04 0.451E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

37 

51 

2 

0 . 132E-02 

0 . OOOE+00 

37 0.462E+04 0.462E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

37 

51 

2 

0 . 150E-02 

0 . OOOE+00 

38 0.410E+04 0.410E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

38 

51 

2 

0 . 189E-02 

0 . OOOE+00 

39 0.417E+04 0.417E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

38 

51 

2 

0 .231E-02 

0 . OOOE+00 


6 . 153947664546494E-3, -1 . 176517428103984E-2 , 2 . 120075657348508E-2 

right left front back 

0.188E+00 0.176E+00 0.223E+00 0.202E+00 

10 trim iterations 

AO =5.946114109135635 A1 =-1.965499082583975 B1 =1.526974903672006 


40 0.339E+04 0.339E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

38 

51 

2 

0 .288E-02 

0 . OOOE+00 

41 0.372E+04 0.372E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

39 

51 

2 

0.369E-02 

0 . OOOE+00 

42 0.401E+04 0.401E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

39 

51 

2 

0 . 515E-02 

0 . OOOE+00 

43 0.404E+04 0.404E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

39 

51 

2 

0 . 657E-02 

0 . OOOE+00 

44 0.396E+04 0.396E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

39 

51 

2 

0 . 848E-02 

0 . OOOE+00 


6.360627238996691E-3, -8 . 764058811065567E-3, 1 . 653152708428962E-2 

right left front back 

0.193E+00 0.184E+00 0.228E+00 0.211E+00 

8 trim iterations 

AO =5.983084257049086 A1 =-2.120740322853678 B1 =1.61160819233671 

45 0.327E+04 0.327E+02 29 39 51 2 0.105E-01 O.OOOE+00 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

{After 45 iterations, the trim routines are able to trim the rotor in less than the 10 iteration limit 
set at the input. The values for Cx and hub moment shown above were ealeulated before the trim 
routines were ealled. } 


46 0.228E+04 0.228E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

39 

51 

2 

0 . 128E-01 

0 . OOOE+00 

47 0.200E+04 0.200E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

40 

51 

2 

0 . 157E-01 

0 . OOOE+00 

48 0.262E+04 0.262E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

40 

51 

2 

0 . 186E-01 

0 . OOOE+00 

49 0.335E+04 0.335E+02 29 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

40 

51 

2 

0 .214E-01 

0 . OOOE+00 


6.458302075405431E-3, -6 . 877644627512502E-3, 1 . 158253065751679E-2 

right left front back 

0.195E+00 0.188E+00 0.228E+00 0.217E+00 

8 trim iterations 

AO =5.969781957344964 A1 =-2.226282329385924 B1 =1.66012949647596 

50 0.394E+04 0.394E+02 29 40 51 2 0.238E-01 O.OOOE+00 

O.OOOE+00 O.OOOE+00 O.OOOE+00 

■k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 

<<<<Restart file has been written>>>> 

(The restart file /ort.20 also eontains the information needed to restart the rotor solution. There 
are no speeial flags that need to be set in the rotor input to handle a restart. } 

CPU time: Total= 2.241E+03; sec/iter= 4.439E+01; sec/iter/point= 108. 7E- 
06 
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